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Chemical structure has been long recognized to greatly influence polymer glass formation, but a 
general molecular theory that predicts how chemical structure determines the properties of glass¬ 
forming polymers has been slow to develop. While the generalized entropy theory (GET) explains 
the influence of various molecular details on polymer glass formation, the application of the GET 
has heretofore been limited to the use of the simplest polymer model in which all united atom groups 
within the monomers of a species interact with a common monomer averaged van der Waals energy. 

However, energetic heterogeneities are ubiquitous within the monomers of real polymers, and their 
implications for polymer glass formation remain to be investigated theoretically. This paper uses an 
extension of the GET to explore the influence of energetic heterogeneities within monomers upon 
the nature of polymer glass formation. This extension of the GET is achieved by combining the 
Adam-Gibbs theory relating the structural relaxation time to the configurational entropy with a 
recent significant extension of the lattice cluster theory for polymer melts with specific interactions, 
in particular, for melts where three distinct van der Waals interaction energies are required to 
describe the energetic heterogeneities within monomers. The present paper focuses on establishing 
general trends for the variation of characteristic properties of glass formation, such as the isobaric 
fragility parameter mp and the glass transition temperature Tg, with molecular details, such as the 
specific interactions and chain stiffness. Our computations confirm that the previously used model 
with monomer averaged interactions correctly captures general trends in the variation of mp and Tg 
with various molecular parameters. More importantly, adjustment of the energetic heterogeneities 
within monomers alone are shown to provide an efficient mechanism for tailoring the properties of 
glass-forming polymers. The variations of polymer properties along iso-fragility and iso-Tg lines are 
illustrated as important design tools for exhibiting the combined influence of specific interactions 
and chain stiffness. 


I. INTRODUCTION 

Despite the fact that polymeric materials readily form 
glasses upon cooling, a general theory for the nature of 
polymer glass formation and for the properties of poly¬ 
mer glasses remains elusive. Numerous studies demon¬ 
strate that molecular characteristics, such as chain stiff¬ 
ness, monomer structure, and the chemical structures of 
the backbone and side groups, profoundly affect poly¬ 
mer glass formation Hence, understanding of the 

relation between material properties and these molecu¬ 
lar details offers the potential to control polymer glass 
formation in a systematic manner. In particular, recent 
experiments 0 indicate that both the glass transition 
temperature Tg and the fragility parameter to, where m 
measures the sensitivity of the structural relaxation time 
or viscosity to temperature changes, can be greatly tuned 
by modifying the chemical structure of the backbone and 
side groups and/or by controlling the spatial positions of 
the side groups with respect to the backbone and/or each 
other. Prior studies illustrate the important changes in 
properties that accompany variations in monomer struc¬ 
ture and thus emphasize the importance of using knowl¬ 
edge concerning the dependence of polymer properties on 
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molecular details to assist in the rational design of poly¬ 
mer materials. However, a complete fundamental under¬ 
standing of how molecular details influence the properties 
of glass-forming polymers remains a challenge. 

The generalized entropy theory (GET) [l^ . [l3| merges 
the lattice cluster theory (LCT) for the thermody¬ 
namics of semiflexible polymers with the Adam-Gibbs 
(AG) theory [H [H , which invokes a relation between 
the structural relaxation time and the configurational 
entropy. Since the LGT employs an extended lattice 
model with structured monomers and describes the in¬ 
fluence on polymer properties of short-range correlations 
imparted by chain connectivity, semiflexibility and inter¬ 
actions, the GET provides a convenient vehicle for sys¬ 
tematically studying the changes introduced in polymer 
glass formation by the alterations of various molecular 
factors. While the agreement of the GET predictions 
with experiment for several nontrivial systems provides 
strong validation of the theory, the goal of rational design 
of polymeric materials requires considering the additional 
complexities of real polymer materials. In particular, all 
previous calculations within the GET [I^ [3 as¬ 

sume all united atom groups within an individual poly¬ 
mer species to interact with the same monomer averaged 
interaction energy. Although this monomer averaged in¬ 
teraction model suffices in establishing general trends 
observed in real polymers, real specific polymers gener¬ 
ally contain energetic heterogeneities within monomers, 
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i.e., different groups have disparate and specific inter¬ 
action strengths. The implication of the influence on 
polymer properties and glass formation of energetic het¬ 
erogeneities within monomers remains to be investigated 
within the GET. 

Motivated by some recent experimental results of 
Sokolov and co-workers 0 and by a desire to minimize 
the enormous complexity of the requisite LCT compu¬ 
tations, we recently extended the LCT to treat a model 
of semiflexible polymers with the structures of poly(n-Q;- 
olefin) and with interactions where only the united atom 
groups at the end of side chains are assigned different 
nearest neighbor van der Waals interaction energies [ 1 ^ . 
Thus, three interaction energy parameters are used to 
describe the specific interactions in this new model. The 
focus of the present paper is to explore how the variation 
of the specific interactions can be used to exert greater 
control over the properties of glass-forming polymers. To 
this end, the new extension of the LCT is combined with 
the AG theory [Slil, thereby providing a similar gener¬ 
alization of the GET to describe polymer glass formation. 
Moreover, the greater realism introduced into the LCT 
and the GET by the new physical model enables testing 
the limits of validity of the simpler LCT model with a sin- 

S monomer averaged van der Waals energy BBS- 
, which is also discussed in detail in the present paper. 


II. POLYMER MELTS WITH SPECIFIC 
INTERACTIONS AND GENERALIZED 
ENTROPY THEORY 


This section presents a brief introduction to necessary 
background information concerning the model of polymer 
melts with specific interactions and the GET of polymer 
glass formation. 


A. Model Polymer Melts with Specific Interactions 

Our new lattice model considers polymer melts where 
the chains are endowed with the structures of poly(n-Q;- 
olefin) (Figure 1 ). The united atom groups residing at 
the ends of the side chains (represented as solid circles 
in Figure 1 and called e-groups in the following) differ 
from those lying on other positions (depicted by open 
circles in Figure 1 and termed n-groups in the follow¬ 
ing), since these side chains’ end segments have different 
attractive nearest neighbor van der Waals interaction en¬ 
ergies among themselves and with other n-groups. The 
microscopic van der Waals interaction energies en and 
£22 are assigned to the interaction of two nearest neighbor 
n-groups and two nearest neighbor e-groups, respectively, 
and the energy ei2 describes the interaction strength be¬ 
tween a n-group and an e-group. This molecular struc¬ 
ture is designed as a minimal model to account for the en¬ 
ergetic heterogeneities within monomers of real polymers. 
These systems are thus termed polymer melts with spe- 
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FIG. 1. Lattice model considered here for the polymer chains 
in melts with specific interactions. Circles designate united 
atom groups, while lines represent the bonds between united 
atom groups. The chains are endowed with the structures of 
poly(n-a-olefin). The example depicted contains side groups 
each containing n = 3 united atom groups. The united atom 
groups residing at the ends of each side chain have specific 
interactions, i.e., different nearest neighbor van der Waals in¬ 
teraction energies £22 with each other and ei 2 with all other 
united atom groups. Therefore, the model prescribes three 
different van der Waals interaction energy parameters (en, 
£22 and £ 12 ). 


cific interactions to distinguish them from polymer melts 
with monomer averaged interactions, where a single van 
der Waals energy specihes the interactions of the system. 
The model for polymer melts with specihc interactions, of 
course, reduces to that with monomer averaged interac¬ 
tions when all nearest neighbor van der Waals interaction 
energies become identical, i.e., en = £22 = £12- 

We view this model of polymer chains with non- 
uniform intermolecular interactions as an idealization of 
polar polymers in which dipolar (or hydrogen bonding) 
species lie pendant off the chain backbone of the polymer 
chain and in which the remaining atomic species along 
the chain backbone interact with only van der Waals in¬ 
teractions. The range of dipole interactions should be 
limited in a medium having a relatively large dielectric 
constant, so that this model can be expected to capture 
some thermodynamic features of these polymer materials 
which have been of great interest with regards to appli¬ 
cations for energy storage devices @. This model also 
serves as a minimal theoretical model for exploring the 
competition between the interactions between two atomic 
species on the nature of glass formation in the GET 
model, a matter of fundamental interest in understanding 
the nature of glass formation since real molecules ordi¬ 
narily have such frustration in interactions. 

In addition to the specific interactions, the LGT also 
includes a description of chain stiffness which is repre¬ 
sented by the presence of a bending energy Ef,. By con¬ 
struction, chains are fully flexible for Eh = 0 , while they 
become completely rigid in the limit Ab —> 00. Also note 
that the side chains may have a different bending energy 
Eg when the number of the united atom groups in each 
side chain n exceeds or equals to two [^. Therefore, 
the free energy derived for the new model is a function 
of polymer volume fraction cj), temperature T, interac¬ 
tion energies (cn, £22, £12), bending energies {Et, Eg), 
molecular weight M, and a set of geometrical indices 
that reflect the size, shape and bonding patterns of the 
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monomers. The technical details of the theory and the 
explicit form for the LCT Helmholtz free energy are pre¬ 
sented in ref for the model of polymer melts with 
specific interactions. 


B. Generalized Entropy Theory of Polymer Glass 
Formation 

The GET treats polymer glass formation as a broad 
transition with four characteristic temperatures M- 
These characteristic temperatures are evaluated from the 
LCT configurational entropy density Sc (i.e., the con¬ 
figurational entropy per lattice site) at constant pres¬ 
sure. This Sc exhibits a maximum as a function of T 
at constant pressure, an essential feature for use in the 
AG model. Recent computations also indicate that 
the LCT configurational entropy density Sc derived in 
ref is nearly identical to the ordinary entropy den¬ 
sity s = —df/dTld, (with / designating the specific 
Helmholtz free energy) for the same thermodynamic con¬ 
ditions, probably because the lattice model is essentially 
devoid of vibrational contributions. Since the ordinary 
entropy density very closely approximates the configura¬ 
tional entropy density and since the former is much 
easier to calculate within the LCT, all calculations em¬ 
ploy the ordinary entropy density in the present paper. 
For simplicity, the ordinary entropy density is just called 
the entropy density in the following. 

The LCT computations for the temperature depen¬ 
dence of the entropy density s(T) enable the direct de¬ 
termination of three characteristic temperatures of glass 
formation, namely, the “ideal” glass transition tempera¬ 
ture To where s extrapolates to zero, the onset temper¬ 
ature Ta which signals the onset of non-Arrhenius be¬ 
havior of the relaxation time and which is found from 
the maximum in s(T), and the crossover temperature 
T/ which separates two temperature regimes with qual¬ 
itatively different dependences of the relaxation time on 
temperature and which is evaluated from the inflection 
point in Ts{T). The conventional definition of the fourth 
characteristic temperature, i.e., the glass transition tem¬ 
perature Tg, requires knowledge of the temperature de¬ 
pendence of the relaxation time r. For this purpose, the 
GET invokes the AG relation ©[il, 

r = rooexp[/ 3 A/rs*/s(r)], ( 1 ) 

where Too is the high temperature limit of the relaxation 
time, /3 = l/iksT) with being ks Boltzmann’s constant, 
A/i is the limiting temperature independent activation 
energy at high temperatures, and s* is the high temper¬ 
ature limit of s(T). Too is set to be 10 “^^ s in the GET, 
which is a typical value for polymers [ 2 ^. Motivated by 
the experimental data on the crossover temperature of 
various glass formers [ 2 ^, the GET estimates the high 
temperature activation energy from the empirical rela¬ 
tion A/i = GksTj [l3l |. Thus, the relaxation time is com¬ 
puted within the GET without adjustable parameters be¬ 


yond those used in the LGT for the thermodynamics of 
semiflexible polymers. The GET then identifies Tg using 
the common empirical definition riTg) = 100 s. In line 
with the original AG theory [I^, the relaxation time r 
calculated from the GET reflects the slowest segmental 
relaxation process (i.e., the primary a-relaxation process) 
in polymer fluids. Other relaxation processes, such as the 
local / 3 -relaxation [30l| . become the dominant relaxation 
below Tg, but they cannot yet be addressed in the GET. 

Once the temperature dependence of the relaxation 
time is known, other related quantities are readily eval¬ 
uated from the GET. For instance, the isobaric fragility 
parameter mp is determined from the standard defini¬ 
tion [m, 


mp 


diogir) 

diTg/T) 


P,T=Tg 


( 2 ) 


Illustrative computations of characteristic tempera¬ 
tures and fragility parameter appear in refs [H, [ 13 , ■ 


III. SOME ILLUSTRATIVE COMPUTATIONS 

This section first discusses the choices for the hetero¬ 
contact (or mixed) interaction energy ei2 for the model of 
polymer melts with specihc interactions. The influence of 
the side group length on glass formation is then explored 
in order to demonstrate the physical validity of the LGT 
model for polymer melts with specific interactions. Fi¬ 
nally, a correlation between the fragility parameter and 
the ratio of different characteristic temperatures, first re¬ 
vealed for the monomer averaged interaction model in 
ref [2^, is shown to likewise apply for the specific inter¬ 
action model. 


A. Influence of the Hetero-Contact Interaction 
Energy 

As discussed in Subsection H A, the model of poly¬ 
mer melts with specific interactions contains three inter¬ 
action energy parameters (en, £22 and £12). Although 
the LGT allows the calculations to be performed with in¬ 
dependent variations of each parameter, analyses of ex¬ 
perimental data suggest that the hetero-contact interac¬ 
tion energy £12 is subject to some constraints and cannot 
be completely independent of £11 and £22. For example, 
it is obviously impossible to hnd a real polymer where 
the hetero-contact interaction energy £12 is significantly 
smaller or larger than either en or £22. Moreover, in¬ 
spired by studies for polymer blends [ 32 |,[ 33 ], an exchange 
energy £ex can be similarly defined for the specific inter¬ 
action model, 


Cex — Cll + £22 — 2£i2. (3) 

The physical essence of £ex for melts with specific inter¬ 
actions is of course different from that for blends. For 
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FIG. 2. Glass transition temperature Tg and isobaric fragility 
parameter mp as a function of 622 for three values of g. The 
computations assume a constant pressure of P = 1 atm, and 
the chains possess the structure of poly (propylene) (PP) with 
molecular weight M = 24001, en = 200 K, and bending 
energy Eh = 400 K. The lines are a guide to the eye. The 
inset depicts the lattice model of polymer melts with specific 
interactions and with the structure of PP, a model which is 
extensively used in the present work. 

instance, the sign of Cex is crucial for determining the na¬ 
ture of the phase behavior of a polymer blend [s^. A 
positive exchange energy implies that the blend phase 
separates upon cooling. However, phase separation ob¬ 
viously cannot occur in a polymer melt. Nevertheless, it 
is instructive to explore whether the sign of Cex affects 
glass formation in the specific interaction model of poly¬ 
mer melts. 

The above considerations suggest two different choices 
for considering the influence of the hetero-contact inter¬ 
action energy ei 2 , i.e., 

£12 = V eiie22, (4) 

and 

£12 = £11 + £22 — •\/£ii£ 22 - (5) 

Both choices limit the number of free interaction energy 
parameters from three to two. The former rule (i.e., the 
Lorentz-Berthelot geometric mean approximation) guar¬ 
antees that the exchange energy Cex is non-negative, while 
the latter sets Cex to be non-positive. We find that both 
choices produce very similar results of Tg and mp (data 
not shown). 

The recent work of Lipson and co-workers [s^ 
demonstrates that a modified geometric mean approxi¬ 
mation for £ 12 , 

£12 = ffV £ii£ 22 , (6) 

provides a better description for understanding the phase 
behavior of polymer blends, where the parameter g char¬ 
acterizes the departure of £12 with respect to the ge¬ 
ometric mean. For instance, a polymer blend is sug¬ 
gested to exhibit a lower critical solution temperature 


for g > 1, while an upper critical solution temperature 
is more likely to occur for g < 1 [H, [s^. Therefore, 
it is interesting to explore the influence of £12 given by 
eq 6 on glass formation in polymer melts with specific 
interactions. 

We consider polymer melts of chains with the struc¬ 
ture of poly(propylene) (PP, where the side group length 
is n = 1) where a single bending energy for the backbone 
Et describes the flexibility of the chains. All calculations 
in the present paper use the same parameter set: the lat¬ 
tice coordination number is 2 ; = 6; the cell volume param¬ 
eter is Oceii = 2.7A; the pressure is P = 1 atm; and the 
molecular weight (i.e., the total number of united atom 
groups in a single chain) is chosen to be M = 24001, cor¬ 
responding to a polymer melt of chains with high molecu¬ 
lar weight. Figure 2 displays Tg and mp as a function £22 
for g = 0.95, 1 and 1.05 [^, when £11 is fixed to be 200 
K. An increase in g leads to an elevated Tg and a slightly 
diminished mp. Nevertheless, the qualitative trends of 
both Tg and mp are not affected by g. Thus, changing 
the sign of (g — 1) has merely quantitative effects on both 
Tg and m p for the specific interaction model. Therefore, 
the subsequent calculations employ the geometric mean 
approximation (i.e., eq 4) to fix the value of £ 12 , with 
the recognition that other choices would yield similar re¬ 
sults. The influence of the interaction energy parameters 
on both Tg and mp is analyzed in detail in Secs. IV and 
V. 


B. Influence of the Side Group Length 

Previous studies demonstrate that the side group 
length significantly influences the glass-formation prop¬ 
erties in polymer melts with monomer averaged inter¬ 
actions [H, [3, m, [31 1 implying that controlling side 
group structure provides a powerful means to regulate 
the properties of glass-forming polymers. Figure 3 dis¬ 
plays both Tg and mp as strongly depending on the side 
group length n in polymer melts with specific interac¬ 
tions. The influence of n on Tg and mp is even stronger 
than that of interaction energy parameters; e.g., the mag¬ 
nitude of Tg at £22 = 100 K for n = 5 is larger than that 
at £22 = 300 K for n = 1. Another noticeable feature in 
Figure 3 is that the change of Tg and mp with £22 weak¬ 
ens for larger n. This result is not surprising since the 
fraction of the side chains’ end segments in a single chain 
varies with n roughly as fe = l/(n -|- 2), indicating that 
the influence of £22 is less important for larger n. Thus, 
Figure 3 also provides evidence for the physical validity 
of the LCT free energy derived for polymer melts with 
specific interactions. 

Notice that Tg increases as the side group length n 
grows for fixed £22 in Figure 3a, a trend that arises be¬ 
cause the illustrative calculations in Figure 3 employ the 
flexible-flexible (F-F) polymer model [31 j where both 
chain backbone and side groups are flexible {Ef, = Eg = 
400 K). A trend of increasing Tg with n occurs also for 
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FIG. 3. (a) Glass transition temperature Tg and (b) iso- 

baric fragility parameter mp as a function of £22 for various 
side group lengths n. The computations are performed for 
polymer melts with specific interactions at a constant pres¬ 
sure of P = 1 atm, where the chain molecular weight is 
M = 24001, the other interaction energies are en = 200 K 
and £12 = ^£ii£ 22 , and the bending energies for backbone 
and side chains are Eh = Ea = 400 K. 


flexible-stiff (F-S) polymers with Et < Eg [I^, a trend 
in accord with measurements for the pair polystyrene 
(PS) 0 and poly (2-vinyl naphthalene) (P2VN) (^, sys¬ 
tems with fairly rigid and extended side groups. Specifi¬ 
cally, experiments indicate that Tg for P2VN exceeds that 
for PS (Tg = 372 K for PS) by 50 K. The GET predicts 
that the other class of polymers, i.e., stiff-flexible (S-F) 
polymers with Eb > Eg, behave in the opposite fashion 
where elevating n leads to a diminished Tg, [13, a 
trend that has been found in polymers with stiff back¬ 
bone and flexible side groups; e.g., Tg is shown in ref 
to decrease significantly with increasing side group length 
for the homologous series of poly(n-alkyl methacrylates), 
polymers with chemical structures that are typical of S-F 
polymers. Therefore, the GET predicted trends for the 
n-dependence of Tg are generally in agreement with the 
observations in real polymers. 
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FIG. 4. Gorrelation between the isobaric fragility param¬ 
eter mp and the characteristic temperature ratio ToITa- 
The computations consider polymer melts with specific in¬ 
teractions at a constant pressure of P = 1 atm, where the 
chains possess the structure of PP with the molecular weight 
M = 24001. The range for each parameter is indicated in the 
legend. 


C. Correlation between Fragility Parameter and 
Breadth of Glass Formation 

A previous work [^ . invoking the monomer averaged 
interaction model, demonstrates that a master curve ex¬ 
ists between the isobaric fragility parameter mp and 
the characteristic temperature ratios such as Tq/Ta- 
This result indicates that the commonly used mp in¬ 
deed correlates with ratios of the characteristic tempera¬ 
tures, supporting the contention that the breadth of the 
glass-formation process provides a promising measure for 
the fragility of glass-forming liquids, [l^ Here, we test 
whether the specific interaction model exhibits similar 
features. Figure 4 displays the correlation between mp 
and To/Ta for independent variations of each parame¬ 
ter (eii, € 22 , or Eb). The collapse of all data, includ¬ 
ing those from the monomer averaged interaction model 
(i.e., squares in Figure 4), indicates that the same master 
curve applies equally for the monomer averaged interac¬ 
tion model [m and for the model of polymer melts with 
specific interactions. Therefore, Figure 4 suggests that 
the fragility generally correlates with the breadth of glass 
formation in glass-forming polymers. 

Next, we focus on how molecular factors control the 
properties of glass formation of melts with specific inter¬ 
actions. Specifically, sections 4 and 5 demonstrate that 
both mp and Tg can be finely tailored by altering the 
molecular parameters, such as interaction energies and 
chain stiffness. Since the simplest model for poly(n-a- 
olefins) with n > 1 contain separate bending energies 
for the backbone and side groups, to simplify the dis¬ 
cussion, the following calculations only consider chains 
with the structure of PP because this choice requires 
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the minimal number of parameters in the LCT [^ . 
Given the chain structure as that of PP, the micro¬ 
scopic interaction energies en, £ 22 , and £12 precisely cor¬ 
respond to backbone-backbone, side group-side group, 
and backbone-side group interaction energies. 


IV. CONTROLLING POLYMER GLASS 
FORMATION BY ALTERING SIDE 

GROUP-SIDE GROUP INTERACTIONS AND 
CHAIN STIFFNESS 

This section illustrates how the properties of glass¬ 
forming polymer melts can be systematically controlled 
by adjusting side group-side group interactions and chain 
stiffness. Since a recent work investigates the com¬ 
bined influence of similar molecular parameters for the 
model of polymer melts with monomer averaged inter¬ 
actions, we test whether the qualitative features found 
for the monomer averaged interaction model apply for 
the more realistic model with specific interactions. This 
section also discusses the behavior of properties along 
the iso-fragility and iso-Tg lines in the plane of the side 
group-side group interaction energy and bending energy, 
concepts first introduced in ref |25| . 

The backbone-backbone interaction energy remains at 
£11 = 200 K for the next set of calculations and hence, the 
side group-side group interaction energy £22 and bend¬ 
ing energy Ei, are the control parameters. This class 
of models can be regarded as polymers with similar 
chemical backbone structures but different side groups 
and is motivated by the recent experimental results of 
Sokolov and co-workers 0 , who explore the alteration of 
polymer glass formation with changes in polar interac¬ 
tions by comparing the properties for polymers with the 
same backbone but different side groups. For example, 
PP, poly(vinyl alcohol) (PVA), and poly (vinyl chloride) 
(PVC) have the same backbone but different side groups. 
Comparing their behavior of glass formation leads to a 
better understanding of how chemical structure controls 
the properties of glass-forming polymers 0. 


A. Combined Influence of Side Group-Side Group 
Interactions and Chain Stiffness 

Figure 5 displays the isobaric fragility parameter mp 
and the glass transition temperature Tg as a function of 
bending energy Ei, for various side group-side group inter¬ 
action energies £ 22 . Both mp and Tg increase with Ef, and 
tend to saturate for very large Ei, at fixed £22 ■ The trends 
of elevating mp and Tg with Ef, are generally in accord 
with physical intuition and experimental observations in¬ 
dicating that greater chain rigidity leads to larger fragili¬ 
ties and higher Tg. For example, mp and Tg for a stiff 
polymer, such as poly(t-butylstyrene) {mp = 141 and 
Tg = 407 K), @ can greatly exceed those for a more flex¬ 
ible polymer, such as poly(dimethylsiloxane) {mp = 85 
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FIG. 5. (a) Isobaric fragility parameter mp and (b) glass 

transition temperature Tg as a function of the bending energy 
Eb for various side group-side group interaction energies £ 22 . 
Solid lines (a) and dotted lines (b) are fits according to eqs 
7 and 8 with the parameters oq = 82.5421, ai = —0.740738, 
02 = 1.10301 X 10■^ bo = 0.412431, for = -1.45264 x 10“^, 
Co = 1.11504 X 10"®, Cl = 2.34999 x 10"'^ and wo = 88.4981, 
ui = -1.07334 X lO'^, U 2 = -3.25406 x 10®, vo = 0.556738, 
ui = 30.3834, V 2 = 3.30207 x 10®, wq = 3.28156 x 10-'‘, 
wi = 0.16638, W2 = 8.24995. 


and Tg = 143 K) [^. Increasing £ 22 , however, leads to 
a drop in mp but an elevation in Tg at the same Eb. 
These qualitative trends are the same as observed for 
the monomer averaged interaction model [l^, [ 2 ^, pro¬ 
viding evidence that the monomer averaged interaction 
model correctly captures general trends in the variation 
of both quantities with interaction energy and chain stiff¬ 
ness. Moreover, two simple algebraic equations fairly ac¬ 
curately capture the computed combined variations of 
mp and Tg with £22 and Eb, 


mp 


0-0 + 0x622 + O2622 + {bo + bx622)Eb 
1 -I- (co + Cl622)Eb 


(7) 


_ Up + m/£22 + U2l^2 + (^0 + ^l/e22 + V2!^2)Eh 

® 1 + (wo + Wl/e22 + W2/ei2)-®b 

( 8 ) 

where the fitted parameters Oa{a = 0,..., 2), ba{a = 0,1), 
Ca{a = 0,1), Ua{a = 0, ...,2), Va{a = 0, ...,2) and 
Wa{a = 0,..., 2) are summarized in the caption of Figure 
5. The forms of the equations are essentially the same as 
those proposed for analyzing the influence of monomer 
averaged interaction e and bending Eb energies on mp 
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FIG. 6. (a) Iso-fragility lines and (b) iso-Tg lines in the plane 
of side group-side group interaction energy 622 and bending 
energy Et for several representative values of mp and Tg. 
Solid lines (a) and dotted lines (b) are the results of eqs 7 and 
8 with the htted parameters given in the caption of Figure 5. 


and Tg for the model of melts with monomer averaged 
interactions [ 2 ^. The fitting functions are chosen to en¬ 
sure the observed saturation of mp and Tg for large Ei,. 

Following the previous work [^ . we define two types 
of special lines in the e22-7?& plane, along which either the 
fragility parameter mp or the glass transition tempera¬ 
ture Tg remains constant. Their existence is due to the 
combined influence of €22 and Eb on mp and Tg. These 
lines are termed iso-fragility and iso-Tg lines, respectively, 
and have been analyzed extensively for the monomer av¬ 
eraged interaction model in order to better understand 
some experimental results of Sokolov and co-workers Q . 

Figure 6 displays several iso-fragility and iso-Tg lines 
as curves of £22 vs. Eb for representative values oi mp 
and Tg, respectively. The lines in Figure 6 display the 
fits from eqs 7 and 8 with the parameters given in the 
caption of Figure 5. In line with the trends observed for 
the monomer averaged interaction model [ 2 ^, Figure 6a 
indicates that Eb grows approximately linearly with £22 
along an iso-fragility line and that the slope grows with 
mp, and Figure 6b exhibits the trend that Eb decreases 
with £22 along the iso-Tg lines, so chains must be more 
flexible at larger cohesive energies in order for the system 
to achieve the same Tg. The trends in Figure 6 can be 
understood as dictated by the combined influence of £22 


and Eb on mp and Tg. The properties along the iso¬ 
fragility and iso-Tg lines are analyzed in subsection 4.2. 


B. Properties along Iso-Fragility and Iso-Tg Lines 

Our previous analysis [ 2 ^ for the monomer averaged 
interaction model demonstrates that many properties, 
e.g., the entropy density, the polymer volume fraction 
and the relaxation time at characteristic temperatures, 
such as Tg, remain invariant along the iso-fragility lines 
in the e-Eb plane. By contrast, no special characteristics 
are found along the iso-Tg lines in the e-Eb plane. There¬ 
fore, the concept of fragility appears to provide more fun¬ 
damental insight into glass formation than the glass tran¬ 
sition temperature Tg. Here, we test whether the model 
of polymer melts with specific interactions displays the 
same properties. 

Figure 7a displays the T-dependence of the entropy 
density s/kp for different pairs of cohesive energies £22 
and bending energies Eb that lie along a typical iso- 
fragility line for mp = 100. The magnitudes of the en¬ 
tropy density at each characteristic temperature remain 
approximately constant along the iso-fragility line in the 
e 22 -Eb plane. This result mirrors that observed for the 
monomer averaged interaction model along iso-fragility 
lines when using the monomer averaged interaction e and 
bending energies Eb as control parameters [ 2 ^, thereby 
providing another evidence for the general utility of the 
monomer averaged interaction model. 

The above result suggests that the entropy density 
along the iso-fragility lines may be a unique function of 
Tq,/T, i.e., the inverse temperature 1/T scaled by one 
of the four characteristic temperatures Tq,. This is con¬ 
firmed by Figure 7b, where the entropy density along the 
same iso-fragility line is presented in an Angell plot. A 
master function approximately describes the dependence 
on Tg/T of the entropy density along the iso-fragility 
lines, although the three curves in the main plot of Figure 
7b slightly deviate from each other at high temperatures. 
The upper inset to Figure 7b also indicates that the re¬ 
laxation times collapse onto a single curve when they are 
plotted as a function of Tg/T, a result that is just a con¬ 
sequence of the iso-fragility lines. By contrast, Figure 8 
reveals that the temperature dependence of the entropy 
density behaves in a rather complicated fashion along 
the iso-Tg lines. In fact, no characteristic properties are 
detected along the iso-Tg lines. The inset to Figure 8b 
displays the monotonic variation of mp with £22 along 
the iso-Tg lines. This behavior accords with many ex¬ 
perimental data, implying that polymers with similar Tg 
may have very different fragilities.For example, PVC and 
poly(3-chlorostyrene) (P3C1S) display similar glass tran¬ 
sition temperatures (Tg = 352 K for PVC and Tg = 362 
K for P3C1S) Q, but their fragilities are quite different 
{mp = 191 for PVC vs. mp = 85 for P3C1S) 0. 

The lower inset to Figure 7b exhibits the glass transi¬ 
tion temperature Tg as increasing with £22 along the iso- 
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FIG. 7. (a) Entropy density s/ks as a function of T for several 
pairs of 622 and Eb that produce the same isobaric fragility 
parameter of mp = 100. (b) Tg-scaled Arrhenius plot for the 
entropy density s/fcp for the same pairs of £22 and Eb as in 
part a. The upper inset to part b presents the Tg-scaled Ar¬ 
rhenius plot for the relaxation time r for the same pairs of £22 
and Eb as in part a, while the lower inset to part b depicts 
£ 22 -dependence of Tg along an iso-fragility line for mp = 100 . 
Squares, circles, diamonds and pentagons designate the po¬ 
sitions of characteristic temperatures Ta, Tj, Tg and To, re¬ 
spectively. Note that the vanishing of the entropy density 
is probably an artifact of high T expansion in the LCT; see 
refs [13, [ 23 ] for discussion. 


fragility lines, a conclusion that is alrea^ evident from 
Figure 7a. Likewise, experimental data [3] indicate that 
the two polymers, PVA and PVC, exhibit very similar 
fragilities (mp = 190), while Tg differs (Tg = 304 K and 
352 K for PVA and PVC, respectively), illustrating that 
there are polymers with similar Tg but different fragili¬ 
ties. The increase of Tg with £22 arises because Eh in¬ 
creases with £22 along the iso-fragility lines in the e 22 -Eh 
plane (Figure 6a). A simultaneous increase of £22 and Eh 
leads to the elevation of all characteristic temperatures. 
Moreover, the lower inset to Figure 7b further reveals 
that Tg increases linearly with £22 along the iso-fragility 
lines in the e 22 -Eh plane. Similar behavior also holds for 
the other characteristic temperatures (data not shown). 

Our previous work using the monomer averaged 



T 



Tg/T 


FIG. 8 . (a) Entropy density s/kp as a function of T for sev¬ 
eral pairs of £22 and Eb that produce the same glass transition 
temperature of Tg = 250 K. (b) Tg-scaled Arrhenius plot for 
the relaxation time r for the same pairs of £22 and Eb as in 
part a. The inset to part b displays the variation of mp with 
£22 along an iso-Tg line for Tg = 250 K. Squares, circles, dia¬ 
monds and pentagons designate the positions of characteristic 
temperatures Ta, Ti, Tg and To, respectively. 


interaction model also implies that the polymer volume 
fraction (j) along the iso-fragility lines in the e-Eh plane 
becomes a unique function of Tg/T and that the volume 
fraction at each characteristic temperature is indepen¬ 
dent of £ and Eh- These properties fail to apply for the 
model of polymer melts with specific interactions, as ev¬ 
idenced in Figure 9a, where the polymer volume fraction 
at Tg is shown to decrease slightly with £22 along the iso¬ 
fragility lines in the e 22 -Eh plane. In addition. Figure 9b 
displays the polymer volume fraction at Tg as increasing 
with £22 along the iso-Tg lines, a trend that is in agree¬ 
ment with that for the monomer averaged interaction 
model [ 2 ^ and can be explained by the negative corre¬ 
lation between £22 and Eh along the iso-Tg lines (Figure 
6b). 
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FIG. 9. Polymer volume fraction at the glass transition tem¬ 
perature (f>Tg as a function of £22 along selected iso-mp and 
iso-Tg lines in the e 22 -i ?6 plane. 


V. CONTROLLING POLYMER GLASS 

FORMATION BY SPECIFIC INTERACTIONS 

Energetic heterogeneities within monomers are ubiqui¬ 
tous in real polymers. Even CHn groups with n = 0,1, 2, 
or 3 are generally assigned different Lennard-Jones inter¬ 
action parameters in continuum model simulations [S- 
Id^ . However, the influence on glass formation of the 
energetic heterogeneities within monomers remains un¬ 
clear. Because modifying the chemical structure of the 
backbone or side groups usually leads to changes not only 
in interaction energies but also in chain stiffness, the in¬ 
fluence of energetic heterogeneities within monomers on 
polymer glass formation cannot readily be isolated ex¬ 
perimentally. The specific interaction model provides an 
opportunity for addressing such issues. Hence, this sec¬ 
tion presents calculations in which the specific interaction 
energies are variables and the bending energy is held con¬ 
stant at El, = 600 K. We demonstrate that the energetic 
heterogeneities with monomers alone can provide an effi¬ 
cient means for tailoring the properties of glass-forming 
polymers. 


A. Combined Influence of Specific Interactions 

Figure 10 displays the isobaric fragility parameter mp 
and the glass transition temperature Tg as a function 
of side group-side group interaction energy €22 for var¬ 


FIG. 10. (a) Isobaric fragility parameter mp and (b) glass 
transition temperature Tg as a function of the interaction en¬ 
ergy £22 for various interaction energies £ 11 . Solid lines (a) and 
dotted lines (b) are fits according to eqs 9 and 10 with the pa¬ 
rameters ao = 316.48, oi = —0.872966, 02 = 1.37215 x 10“®, 
bo = 0.365033, bi = -3.32456 x 10“®, 62 = 6.63404 x 10“®, 
Co = 6.43496 x 10"^ ci = -3.25912 x 10"®, C 2 = 6.81742 x 
10"® and uo = 48.2178, m = 0.913154, U 2 = -9.17928X lO"'^, 
Vo = 0.883654, vi = 3.21211 x 10"®, V 2 = -4.65713 x 10“®, 
wo = 2.65577 x 10“®, wi = 7.49098 x 10"®, W 2 = -1.79918 x 
10 "®. 


ious backbone-backbone interaction energies en. The 
fragility parameter mp decreases with increasing either 
£11 or £ 22 , while elevating either of the interaction ener¬ 
gies leads to an increase in Tg. Apparently, the separate 
influence of £11 or £22 on mp and Tg is qualitatively the 
same as that of the monomer averaged interaction energy 
£. Likewise, we find that two slightly different algebraic 
equations capture the computed combined variations of 
mp and Tg with cn and £ 22 , 


ao + ai£ii -|- a2£ii + [bo + 6i£ii -I- 62 eii)e 22 
=-;- 7 -ry-^ (9) 

1 + (Cq + Ci£ii + C2£ii)£22 


J, _ Up + Uieii + U2e\i -I- (uq -I- ui£ii -I- r;2en)£22 

® 1 + (wo + Wi£ii-I-r(;2eii)e22 

( 10 ) 

where the fitted parameters aa{a = 0,...,2), ba{a = 

0,..., 2), Ca{a = 0,..., 2), Ua{a = 0,..., 2), Va(a = 0,..., 2) 

and Wa{a = 0,..., 2) are provided in the caption of Figure 
10. Figure 10 demonstrates that the energetic hetero¬ 
geneities with monomers alone provide additional vari¬ 
ables for tailoring the properties of glass-forming poly¬ 
mers. 
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FIG. 11. (a) Iso-fragility lines and (b) iso-Tg lines in the plane 
of interaction energies en and £22 for several representative 
values of mp and Tg. Solid lines (a) and dotted lines (b) 
are the results of eqs 9 and 10, respectively, with the fitted 
parameters given in the caption of Figure 10. 

The iso-fragility and iso-Tg lines can be similarly de¬ 
fined in the £ 11-622 plane. Figure 11 displays several iso¬ 
fragility and iso-Tg lines as curves of en vs. £22 for rep¬ 
resentative values of mp and Tg, respectively. The lines 
in Figure 11 display the fits obtained from eqs 9 and 10 
with the parameters given in the caption of Figure 10. 
The features of these lines in the £ 11-622 plane are clearly 
different from those in the 622 -F'f) plane. This result is ex¬ 
pected since the influence of interaction energies on mp 
and Tg differs from that of bending energy in quantita¬ 
tive or even qualitative ways. Figure 11 indicates that a 
negative correlation exists between £11 and £22 for both 
iso-fragility and iso-Tg lines. Again, such a result can be 
explained by analyzing the combined influence of £11 and 
£22 on mp and Tg. 


B. Properties along Iso-Fragility and Iso-Tg Lines 

We now provide similar analyses for the properties 
along the iso-fragility and iso-Tg lines in the £ 11-622 plane 
and test whether the characteristic features found in the 
622 -F'b plane can be generalized when the bending energy 
is held constant. 

Figure 12a displays the T-dependence of the entropy 
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FIG. 12. (a) Entropy density s/ks as a function of T for 
several pairs of eii and £22 that produce the same isobaric 
fragility parameter of mp = 130. (b) Tg-scaled Arrhenius 

plot for the entropy density s/fcp for the same pairs of £11 
and £22 as in part a. The upper inset to part b presents 
the Tg-scaled Arrhenius plot for the relaxation time r for the 
same pairs of £ii and £22 as in part a, while the lower inset 
to part b depicts £ii-dependence of Tg along an iso-fragility 
line for mp — 130. Squares, circles, diamonds and pentagons 
designate the positions of characteristic temperatures Ta, Tj, 
Tg and To, respectively. 


density for different pairs of interaction energies £11 and 
£22 that lie along an iso-fragility line for mp = 130. In 
contrast to the results in the £ 22 -F'& plane, the magni¬ 
tudes of the entropy density at characteristic tempera¬ 
tures (Ta, Tj, and Tg) clearly increase with growing cn 
or diminishing £22 along the iso-fragility lines in the £ 11 - 
£22 plane. This immediately suggests that the scaling of 
the entropy density breaks down along the iso-fragility 
lines in the £ 11-622 plane (Figure 12b) and implies that 
the features displayed along the iso-fragility lines in the 
£ 22 -A{, plane cannot be fully generalized. Nevertheless, 
the characteristic temperatures still increase linearly with 
£11 along the iso-fragility lines in the £ 11-622 plane (see 
the lower inset to Figure 12b). The upper inset to Fig¬ 
ure 12 b depicts the relaxation times as still collapsing 
onto a single curve when plotted as a function of Tg/T, 
and the relaxation times at the characteristic tempera- 
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FIG. 13. (a) Entropy density s/fcs as a function of T for sev¬ 
eral pairs of en and £22 that produce the same glass transition 
temperature of Tg = 270 K. (b) Tg-scaled Arrhenius plot for 
the relaxation time r for the same pairs of en and £22 as in 
part a. The iirset to part b displays the variation of mp with 
£11 along an iso-Tg line for Tg = 270 K. Squares, circles, dia¬ 
monds and peirtagons designate the positions of characteristic 
temperatures Ta, Ti, Tg and To, respectively. 


ture Ta, Ti, or Tg remain approximately constant along 
the iso-fragility lines in the eii-e 22 plane. Moreover, the 
iso-Tg lines in the eii-e 22 plane are again not found to 
possess special features in the temperature dependence 
of both entropy density and relaxation time (Figure 13). 
The inset to Figure 13b exhibits the monotonic varia¬ 
tion of mp with eii along the iso-Tg lines in the eii-e 22 
plane, which implies that different polymers with similar 
Tg may exhibit large variations in mp even when these 
polymers possess a similar degree of chain stiffness. 


Figure 14a displays the polymer volume fraction at Tg 
as clearly increasing with en along the iso-fragility lines 
in the e-Ei plane and, thus lacking the scaling of the 
volume fraction at Tg obeyed by the simpler monomer 
averaged interaction model 2^. This result emerges be¬ 
cause van der Waals interactions exert a dominant role 
in determining the equation of state (EOS) of a polymer 
melt. [ 2 ^ Hence, altering both the interaction energies 
£11 and £22 is expected to induce complicated changes 



Oi 

FIG. 14. Polymer volume fraction at the glass transition tem¬ 
perature 4>Tg as a function of £11 along selected iso-mp and 
iso-Tg lines in the £ii -£22 plane. 


in the EOS, a result that also applies along the iso-Tg 
lines in Figure 14b, which illustrates the nonmonotonic 
dependence of the polymer volume fraction at Tg on en 
along the iso-Tg lines in the £ii -£22 plane. 


VI. DISCUSSION 

Although energetic heterogeneities within monomers 
are ubiquitous in real polymers, no prior theory can ad¬ 
dress their role in polymer glass formation. The present 
paper extends the GET to consider the role of the ener¬ 
getic heterogeneities with monomers in determining the 
properties of glass-forming polymers. This achievement 
is made possible because of the recent, technically com¬ 
plex extension of the LCT to describe polymer melts with 
specific interactions [ 2 ^ , where three van der Waals inter¬ 
action energy parameters are used to account for the en¬ 
ergetic heterogeneities within monomers in real polymers. 
The greater physical realism introduced into the GET 
enables testing the limits of validity of the monomer av¬ 
eraged interaction model, which has been extensively em¬ 
ployed to investigate the essential and physical molecular 
features affecting polymer glass formation [l^, , 

to model the glass formation of poly(Q;-olefins) , to an¬ 
alyze the thermodynamic scaling of dynamics in polymer 
melts [ 2 ^, and recently to address the two glass tran¬ 
sitions in miscible polymer blends (s^ . [ 3 ^ . The com¬ 
putations confirm that the monomer averaged interac¬ 
tion model correctly captures the general trends of the 
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variation of mp and Tg with various molecular param¬ 
eters. More importantly, the energetic heterogeneities 
with monomers alone are shown to provide an efficient 
way for tailoring the properties of glass-forming poly¬ 
mers. 

Using the extension of the GET, we further explore 
how the characteristic properties of glass formation in 
polymers can be controlled by specific interactions and 
chain stiffness using a minimal model containing a sin¬ 
gle bending energy. We focus on two cases. The first 
case employs a fixed backbone-backbone interaction en¬ 
ergy, and the side group-side group interaction energy 
£22 and bending energy Ei, are used as control param¬ 
eters. This class of models is relevant to understand¬ 
ing the groups of polymers recently studied by Sokolov 
and co-workers 0, who explore the influence of polar 
interactions on polymer glass formation by comparing 
the results for polymers with the same backbone but dif¬ 
ferent polar or hydrogen bonding side groups, the class 
of materials for which our model is an idealization. In 
the second case, the bending energy is held constant and 
hence, the changes in properties are solely caused by the 
specific interactions. Iso-fragility and iso-Tg lines are de¬ 
fined in both cases, and we analyze the properties along 
the iso-fragility and iso-Tg lines defined in both cases. 


Many features displayed by the monomer averaged in¬ 
teraction model also appear along those special lines in 
the model of polymer melts with specific interactions. 
However, some differences are also evident. For exam¬ 
ple, the scaling of the polymer volume fraction observed 
in the monomer averaged interaction model along the 
iso-fragility lines is absent in the model of polymer melts 
with specific interactions. Nevertheless, the results in the 
present paper clearly demonstrate the physical validity 
of the monomer averaged interaction model and provide 
evidence that tailoring molecular details such as specific 
interactions and chain stiffness provides an efficient route 
for guiding the rational design of glassy polymer materi¬ 
als. 
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